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Abstract 

When a color superconductor of high density QCD is rotating, super- 
fluid vortices are inevitably created along the rotation axis. In the 
color-flavor locked phase realized at the asymptotically large chemical 
potential, there appear non-Abelian vortices carrying both circulations 
of superfluid and color magnetic fluxes. A family of solutions has a 
degeneracy characterized by the Nambu-Goldtone modes CP 2 , associ- 
ated with the color-flavor locked symmetry spontaneously broken in the 
vicinity of the vortex. In this paper, we study electromagnetic coupling 
of the non-Abelian vortices and find that the degeneracy is removed 
with the induced effective potential. We obtain one stable vortex solu- 
tion and a family of metastable vortex solutions, both of which carry 
ordinary magnetic fluxes in addition to color magnetic fluxes. We dis- 
cuss quantum mechanical decay of the metastable vortices by quantum 
tunneling, and compare the effective potential with the other known po- 
tentials, the quantum mechanically induced potential and the potential 
induced by the strange quark mass. 



1 Introduction 



The color superconducting phase occupies the low temperature-high density region of the 
QCD phase diagram fil. For these values of physical parameters, quarks condense in couples, 
similarly to what happens in ordinary superconductors for electrons forming Cooper pairs. 
Moreover, this phase is classified with different names, depending on which flavors participate 
to condensation. The first rough distinction is made between the so-called 2SC and color- 
flavor locked (CFL) phase (2j|3]: the first one is characterized by the condensation of the two 
lighter flavors only; in the other all three flavors participate to the condensation. It has been 
discussed that the recently found heavy neutron star |4] implies the equation of state soft, 
which denies the existence of exotic matters including quark matters. However, there is also 
an objection to that argument [H] and the situation is not conclusive yet. The CFL phase 
may exsist in the core of compact stars, and it is important to suggest some experimental 
signatures of its exsistence. 

In the CFL phase, due to the condensation of the three quark flavors, the original sym- 
metry group is broken down as G = U{1)b X SU(3)c x SU(3)f — > H = SU{3)c+f x Z 3 , 
where U(1)b is the baryon symmetry group, while, apart from the color and flavor symme- 
tries, SU(3)c+f is the color-flavor locked group. Since U(1)b and SU(3)c are spontaneously 
broken in the ground state, it is superfluid and color superconducting, respectively, at the 
same time. When a color superconductor is rotating, which will be the case if it is realized in 
a neutron star core, superfluid U(1)b vortices (6]^8j are inevitably created along the rotation 
axis, as in usual superfluids such as 4 He and Bose-Einstein condensates in ultracold atomic 
gases. However a U(l) vortex is unstable to decay into a set of three non-Abelian semi- 
superfluid vortices found by Balachandran, Digal and Matsuura [9 |Q because of a long range 
repulsion between them (12| . Then such non-Abelian vortices are expected to form a vortex 



lattice 12, 13]. The precise numerical solutions and analytic expression of the asymptotic 
form of the vortex solution were given in (l4j. Since a non-Abelian vortex winds 2n/3 of the 
U(1)b phase and SU(3)c at the same time, it carries 1/3 circulation of that of the U(1)b 
vortex and a color magnetic flux, respectively. Non-Abelian vortices are superfluid vortices 
as well as color magnetic flux tubes so that they are the most fundamental vortices in color 
superconductors . 

These solutions are invariant only under the action of a subgroup of the color-flavor 
symmetry of the ground state, thus producing the further breaking SU(3)c+f — > SU(2)c+f x 
U(1)c+f- Consequently, the Nambu-Goldstone modes, called orientational moduli, 

SU(2)o+f x U(l) c+ F ~ 1 ' 



arise from this breaking near the vortex core |12|. There is a whole family of solutions 
all having the same tension which can be mapped the ones into the others by color-flavor 



transformations. It was shown in 15 that the orientational zero modes (1.1) are in fact 
normalizable, and the effective world-sheet theory was constructed as a CP 2 non-linear sigma 
model. The effect of the strange quark mass was taken into account in the effective world- 
sheet theory as the potential term |16|. In the asymptotically large chemical potential where 



1 Non-Abelian vortices were found in the context of supersymmetric theories 10 , in models whose vacuum 



possesses the color- flavor locking structure as the CFL does. See Refs. 11 as a review in this context. 
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strange quark mass can be neglected, quantum effects imply the presence of monopoles, seen 
as kinks on the vortex worldsheet [17] ■ Coupling of non-Abelian vortices to quasi-particles 
such as gluons, U(1)b Nambu-Goldstone bosons (phonons) were determined [18]. More 
recently, the electro-magnetic coupling to the non-Abelian vortices has been studied and it 
has been found that a lattice of non-Abelian vortices acts as a polarizer of photons fi.9|. Some 
implications of the existence of non-Abelian vortices in neutron stars were studied |20) . These 
analysis were all based on the Ginzburg-Landau model which is valid for temperatures near 



the critical one 21,22 . In addition to the analyses based on the Ginzbur-Landau model, the 
Bogoliubov-de Gennes approach which is valid below the critical temperature was also used 
to study the fermionic structures of the non-Abelian vortices, and it was found that non- 
Abelian vortices trap Majorana fermion zero modes belonging to the triplet representation 



of the unbroken symmetry S77(2)c+f at their cores 23 . One of intersting consequences of 
the fermion zero modes is that the exchange statistics of non-Abelian vortices makes them 
non-Abelian anyons |24|. 

In this paper, we study electromagnetic coupling of non-Abelian vortices based on the 
Ginzbur-Landau model. While the analysis in [19] neglected the mixing of the photon and 
the gluon, we correctly taking it into account. Our analysis is, in some sense, a continuation 
of the work done in j9j, where also electromagnetic coupling was considered. Since the Z7(1) BM 
symmetry is embedded into the flavor SU(3)p symmetry, the latter symmetry is not intact 
anymore; only its subgroup SU(2)p x t/(l)p remains exact. We find, as a consequence, that 
it removes the degeneracy of solutions (1.1 ) in energy between the different vortex solutions. 
We show that the solution of |9| (which we call "the BDM vortex") represents the lowest 
energy configuration when we neglect the strange quark mass, which was in fact stated 
in |9|. We are able also to identify two other different "diagonal" vortex solutions having a 
tension slightly higher compared to the former, and being related together by the residual 
SU(2)c+f and then constituting a full family of solutions, a CP 1 space as a submanifold of 
the CP 2 moduli space. We refer to these last solutions as "CP 1 vortices". We find another 
special vortex solution whose winding does not regard the direction in the color space that 
allows for the electromagnetic gauge field to couple, namely the T 8 direction in SU (3)c- This 
solution possesses the highest tension and should be unstable. We call it "pure color" vortex. 
An energy potential connects all the different vortices: the BDM solution is its absolute 
minimum, while the two CP 1 vortices are local minima, along with their full CP 1 space; the 
pure color vortex stays at the maximum. We then find that the BDM and CP 1 vortices are 
spontaneously magnetized to carry magnetic fluxes of the electromagnetic field, in addition 
to the color magnetic fluxes. The magnetic flux of the BDM vortex is —2 times the one of 
the CP 1 vortices. Since the minimal color-less bound state is obtained as the U(1)b vortex 
made of one BDM vortex and two CP 1 vortices, this color-neutral bound state of vortices 
necessarily carries also no electromagnetic flux, and vice versa. The amount of magnetic flux 
is proportional to e 2 /g 2 and is tiny. 

This situation is changed when the strange quark mass is taken into account. In this case, 
as discussed in Ref. [16] , a potential arises between the vortex configurations considered so 
far, breaking the residual CP 1 symmetry and lowering the tension of one of the CP 1 vortices, 
which will be called CP+ vortex. This solution has now the lowest energy and is then the 
most stable one, to which the BDM vortex can decay. We will show that, for densities typical 
of the inner neutron star core, the CP} vortex is the relevant solution, instead of the previous 
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considered BDM vortex. 

This paper is organized as follows. In Section [2] we introduce the theoretical setting in the 
case without electromagnetic coupling and we review the pure color solutions. In Section [3] 
we add the electromagnetic coupling to our model and we explore the different features of 
the three kinds of vortices mentioned above: the BDM vortex, the CP 1 vortices and the 
pure color vortex. We calculate magnetic fluxes which these vortices carry. We then discuss 
the CP 1 vortices, which are metastable classically, decay quantum mechanically into the 
BDM vortex by a quantum tunneling. We give an estimate of the decay probability. Then 
in Section [4] we compare the potential induced by the electromagnetic interaction which 
we found with other potential terms; the one quantum mechanically induced and the one 
induced by a non-zero strange quark mass. In Section [5] we present our conclusions. In 
Appendix, we summarize equations for non-Abelian vortices. 

2 QCD in the CFL phase 

2.1 Symmetry breaking in the CFL phase 

Color-superconductivity in hight density QCD is almost "inevitable". When the chemical 
potential \x is large enough, \i > Aqcd, QCD is in a perturbative regime. In this regime, it 
is easy to see that gauge interactions mediate attractive forces between quarks. The BCS 
mechanism of superconductivity implies then formation of Cooper pairs and of a diquark 
condensate (2J. In the most symmetric phase, (color-flavor-locked, CFL), the diquark con- 
densate has the following form (3]: 

(VfCV/f) = A CPL [e^%- fc $ 7 fc ] , i,j,k = u,d,s a,P,-y = r,g,b, (2.1) 

where the two quarks pair in a parity-even, spin singlet channel, while the color and flavor 
wave functions are completely antisymmetric. The order parameter $ is a 3 by 3 matrix 
that can be written, in terms of pair condensation, as follows 

/ (ds) gb (us) gb (ud) gb \ 
K = (ds) rb (us) rb (ud) rb . (2.2) 

\ (ds) r g ( US )rg ( ud ) rg J 

It transforms under gauge and flavor symmetries as a (3c, 3p) representation : 

$ -> U C $>U ] F . (2.3) 

In the CFL phase, the ground state is given by the following value of the order parameter 
(up to color-flavor rotations): 

<$ 7 fc ) = A CFL 5 7 fc . (2.4) 

We notice en passant that other phases are possible when one considers a non-zero mass 
for the strange quark. For example, when the chemical potential is of order of \i ~ m^/A CFL , 
the strange quark does not participate to condensation anymore and we have the so-called 
2SC phase, where the order parameter is given by $ k = A2 SC 8 1 3 S 3 k . 
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As well-known, QCD has the following symmetries when quark masses are neglected: 



G = U(1) B x SU(3) C x SU(3) L x SU(3) R . (2.5) 



The condensate in Eq. (2.4) breaks them to the following residual group: 

H = SU(3)c+l+r. (2.6) 

The residual non-Abelian symmetry in the ground state is a peculiarity of the CFL ground 
state and the origin of some of the most interesting properties of this phase. Notice for 
example that chiral symmetry is broken perturbatively. 

The symmetry breaking pattern described above implies the existence of stable topolog- 
ical "semi-superfTuid" vortices [9] 

..W^ f"" 1 ^ )^. (2.7) 

The origin of the word semi-superfluid arises because the smallest non-trivial loops corre- 
sponding to vortices involve both global and gauge rotations. Vortices in the CFL phase of 
QCD are similar to vortices appearing in both superfluid and superconductors. 

In the next Section we will briefly review the construction and the properties of these 
vortices. 



2.2 Landau- Ginzburg description of the CFL phase 



The Landau-Ginzburg description of the CFL phase in terms of a field theory for the order 
parameter is appropriate at temperatures close to the critical temperature T c for the CFL 
phase transition. Since we will always be interested in static configurations, we write down 



only the terms that do not include time derivatives 21 , 25 

"eA 



LG 



Tr 



F 0j F 0j - -FijF ij + K 1 V i &V i $ - A 2 ($ t $) 2 + m 

$t$ 



- Ai(Tr [$ T $]) - 



4" 13 

3m 4 
4(3Ai + A 2 ) 



+ 



(2.1 



Traces are taken on both color and flavor indexes when as appropriate. We use the following 
conventions: 



A- 



Fij — d- L Aj — djAi — ig s [Ai, Aj], 
V i <$> = d i <$>-ig s A i $, 
A?T a , F^ = i^-T a , Tr (T a T b ) = 5 ab . 



(2.9) 



The coefficient in the expression above may be calculated directly from the QCD lagrangian 
using perturbative techniques. We quote here the standard results obtained in literature 



5 



through perturbative calculations in QCD [7|[2l]: 



m 2 = -4iV(/i)log-, iV(/i) 



24ir 2 A / 3ir 2 \ , 

9 '='/271^' Tc ~" eXP l"^J' (2 ' 10) 

where /i is the chemical potential, A the QCD scale and T c the critical temperature, while 
e and A are the dielectric and diamagnetic constants. In the following, we assume electric 



fields to be vanishing, and we thus neglect the first term in Eq. (2.8). 



The expectation value of $ can be directly found by minimizing the potential of expres- 



sion (2.8), obtaining: 

,2 



W =A CPL 1 3 , A^^. (2.11! 



Masses of gauge bosons and scalars are given by the following: 

o m 2 a\ A 2 

2 o 2a2 ts 2 2 ^^^cfl 2 n /nm\ 

™ q = 2 9s A cfl K u = -—, m = — , m = 0, (2.12) 



where ip is the massless Nambu-Goldstone boson related to the breaking of U(1)b symmetry, 
and <j) and X are respectively the trace and traceless part of $. 

Equations of motion can also be directly found from the lagrangian above, and they read 

as: 



V, /••'•' = -ig s K 



Vj<S>& - ^(V^) 1 - ^Tr (Vj$$ f - $(Vj$) f ) 



VjV J '$ = J- [-\ 2 $& - 2AiTr ($ f $) + m 2 ] $ . (2.13) 
2.3 "Pure" color-magnetic flux tubes 

Let us discuss now in more detail the color flux tubes that are stabilized by the topology 
of the CFL phase. These solution have been thoroughly studied in literature |9| [l5"j|l2~|[l4] 
and have been dubbed, in various works, as semi-superfluid, color-magnetic or non-Abelian 
strings. In the following we will also refer to these vortices as un- coupled vortices, because 
in Section [3] we will couple them to the electromagnetic gauge field. The latter solutions will 
obviously be the coupled ones. 

Strings in the CFL are semi-superfluid because the fundamental vortex corresponds to 
a non-trivial closed loop constructed with both global U(1)b and gauge SU(3)c symmetry. 
This fact can be easily seen by more closely inspecting the relevant homotopy group in 



Eq. (2.7). The Z3 factor in the denominator implies that the smallest loop can be constructed 
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by picking up a 7r/3 global phase and an additional 27r/3 gauge phase. More concretely, a 
vortex configuration can be written as follows: 



e**7(r) 
g(r) 
g(r) 



$(r,^) f = A CFL | g(r) |, A 8 T 8 = ^ ^[1 - h(r)} T 8 , (2.14) 



where the profile functions /(r), g(r) and /i(r) satisfy the equations reported in the Appendix 
and satisfy the following boundary conditions: 

/(0) = 0, </(0) = 0, h(0) = l, /(oo) = g(oo) = 1, /i(oo)=0. (2.15) 

The vortex above corresponds to the following closed loop in the global-gauge group: 

$(oc, (p) = e^ /3 e" 4 v / 273^ T8 $(oo, 0) . (2.16) 

Since the (ds) gb pairing winds in this case, we called the previous configuration "anti-red" 
vortex. By inserting the asymptotical value above into the original lagrangian (2.8), we can 
extract a logarithmically divergent contribution typical to global vortices. For the configu- 
ration above the tension is: 



T= 2ixK 1 I rrfr Tr I Vj$r ^ 2nK 1 / rdr Tr 
2tt 



A^Ki In L/L core , r oo (2.17) 

where L and L core denote the sizes of the system and the vortex core, respectively. The diver- 
gent contribution to the tension can be put in relation with a generic fractionally quantized 
circulation eg of the superfluid quarks: 

T =—A 2 CFL K 1 In L/ L core = 67rc|A^ FL ifi In Lj L core , 

7TI* f -> 71 

c B =— <t> v sup ■ dl = - . (2.18) 

Semi-superfluid vortices in the CFL phase are also color-magnetic in the sense that they 
also carry a non-trivial color flux. For the configuration above we have: 



Since the vortex configuration (2.14) breaks the SU(3)c+f down to a U(1)c+f> we can 
straightforwardly construct the most general vortex configuration by applying color-flavor 
rotations: 

$(r, ^general = U${r, f^U' 1 , (2.20) 

which gives a full set of degenerate configurations given by the following space: 

M = CP 2 (2 21) 

McFL ~ SU(2) x U(l) ^ ■ 1 ' 
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We can schematically represent the space CP 2 in terms of its "toric" diagram, as in Fig. [l] 
Each point in the internal part of the diagram represents a 2-torus generated by the U(l) 2 
isometry of CP 2 . Each point on the sides represents a torus where a U(l) degenerates. 
The vertices of the triangles represent just points in the moduli space. In our notations, 
the vertices are represented by the 3 diagonal vortex configurations where the winding is 
concentrated on each diagonal entry only, giving respectively f, g and b vortices. 




Figure 1: The "toric diagram" of CP 2 . The different vortices f,g,b are located at the 
vertices: the winding of these vortices is localized in one of the diagonal entries in the 
matrix field <F The upper-right edge of the diagram corresponds to the CP 1 subspace 
invariant under the SU(2)c+f, whose generators commute with T EM defined in Sectionp) 



To conclude this Section, we discuss the dependence of the energy on the coupling con- 
stant g s . In general, the tension of the vortex has a logarithmically divergent contribution 
that depends only on the gap. However, finite corrections to the energy depend in general on 
the gauge coupling and on the various coefficients in the potential energy. We have calculated 
these finite contributions numerically and we have studied their behavior under the change 
of g s . Fig. [2] shows that the tension decreases monotonically with increasing of the gauge 
coupling. The behavior is generic and independent of the values of scalar masses. It can 
be intuitively understood noticing the 1/ 1 g 2 s dependence in the kinetic terms for the gauge 
potentials. 



3 CFL Phase with Electromagnetic Interactions 
3.1 Electromagnetic Interactions 

The introduction of electromagnetic interactions is straightforward once we look at the charge 



structure of the order parameter as in Eq. (2.2). The charges of each column of $ are 



respectively Qem^) = (— 2/3, 1/3, 1/3) and electromagnetic gauge transformations can be 
implemented by a right multiplication generated by: 

$ § e iaTE \ T EM = - diag(-2, 1, 1), T 8 = ^/3/2T EM . (3.1) 

3 
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Figure 2: The figure shows the monotonic decreasing of the corrections to the tension 
of the vortex as function of the gauge coupling, which in this figure is plot against the 
mass of the gauge bosons m g . The red line corresponds to the "realistic" parameters: 
H ~ 500MeV, A ~ 200Mev, T c ~ lOMeV, T ~ 0.9T C , which correspond to A CFL =7 Mev, 
K\ = 9, TO^ = 34Mev, m x = 17Mev. The red and green lines correspond respectively 
to (m^,m x ) = (50Mev,10Mev) and (m^,m x ) = (10Mev,50Mev). The logarithmically 
divergent tension has been cut-off at a distance L = 0.4MeV _1 



We can now generalize the Landau-Ginzburg description of the previous Section: 



£ = Tr 



Spvupmuij _ ^FijF ij + KjVi&V*® - A 2 ($ t $) 2 + m 2 $ f $ 



+ 



where the coefficient in front of the kinetic term for the electromagnetic field is due to the 
different normalization for the T EM generator; we have also used the following conventions: 

Vi$ = di<$> - ig a A@ - ie$Af M . (3.3) 

The ground state of the theory is unchanged with respect to the un-coupled case: 

<$ 7 fc )=A CFL 5 7 fc , (3.4) 

however, the first effect of introducing electromagnetic interactions is the explicit breaking 
of the SU (3) f flavor symmetry, even if we still consider all quarks to be massless. Actually, 
introducing electromagnetic interactions is equivalent to gauging a U(l)p subgroup of the 
SU(3)f flavor symmetry. 

SU(3) F ^SU(2) F xU(l) EM . (3.5) 

The full set of symmetries of the CFL phase of QCD with electromagnetic interactions is 
thus given by the following (we omit chiral symmetry): 

G = U(1) B x U(1) EM x SU(3) C x SU(2) F . (3.6) 

The second important effect of introducing electromagnetic interactions in the CFL 
ground state is the presence of an unbroken U(1)q symmetry which is a combination of 
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T EM and the T 8 generator of SU(3)c- They have to be rewritten in terms of their massless 
and massive combinations [I] 

A M = cos(A em + sm(A 8 ■ 
Ao = -sin(A em + cos(A 8 , (3.7) 

written in terms of a mixing angle ( 



e 2 e 



e 2 + 3^/2 g M 

Notice that the mixing above is really meaningful only when the left action of T 8 and the 
right action of T EM are really equivalent, which means for diagonal configurations of the order 
parameter $. In this case, we can for example write the covariant derivative like: 

Vi^di-iguAf . (3.9) 

Because of this non-trivial mixing, the masses of gluons are not equal in the CFL phase: 

m 2 a = 2g 2 A 2 CFh K u m 2 8 = 2g 2 M A 2 CFJj K 1 , (3.10) 

where a runs from 1 to 7 and mg is the mass of the gluon related to the meted generator T 8 . 



Apart from the unbroken gauge U(l) symmetry, the CFL ground state Eq. (3.4) has the 
following diagonal color-flavor symmetry 

H EM = SU(2) C+ F ■ (3.11) 

This reduced symmetry is crucial to understand the property of the moduli space of non- 
Abelian with electromagnetic coupling. 

3.2 "Coupled" color-magnetic flux tubes 

We approach the study of vortices in the CFL phase when they are coupled to electromagnetic 
fields in the most general way, starting from their topological classification. The most general 
vortex is related to all the possible non-trivial loops in the vacuum manifold. In the CFL, 
this manifold is exclusively generated by symmetry transformations: 

(ka \ <niu\ f U ^ B x SU(3) C -f \ v rqio v 
iri{M va c) =-K 1 {G/H) =7Ti — )~Z. (3.12) 

V [£>3)C-F+B J 

The Z3 factor above is crucial to have semi-superfluid vortices with non-Abelian fluxes. 
This is due to the fact that the smallest non-trivial loop has to wind in both the global 
baryonic and gauge symmetry group. The conclusion above is not changed if we consider 
the additional electromagnetic coupling. This is due to the fact that we can always unwind 
an electromagnetic phase by using the unbroken gauge group C/(l)o- However, in the un- 
coupled case we have a full SU(3) isometry that we can use to generate all the possible vortex 



configurations starting with an arbitrarily chosen one, for example the one in Eq. (2.14). On 
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the other hand, in the coupled case we only have a residual SU(2) isometry, which is not 



enough to exhaust all possible solutions. For example, the vortex in Eq. (2.14) is invariant 
under this symmetry. In the toric representation of Fig. [T] the SU(2) isometry transforms 
points of the triangle along lines parallel to the long diagonal side. See also Fig. [TJ 

Let us make this discussion more concrete. A closed loop in the gauge/baryon group is 
given by the following transformation on the order parameter, where <p is the angle coordinate 
of the space: 

($(oo,y?)> = e ie ^e l ^ ( ^ Ta ($(oo,0)) e ia ^ TEM 

J| (f — 27T 

e i~t a (2-K)T a _ e -i0(27r) e -iQ(27r)T EM ^ (3.13) 

where 6, 7" and a are monototically increasing functions. The second line is an equation for 
the possible symmetry transformations giving a closed loop, or, equivalently, a possible vortex 
configuration. Notice that existence and stability of vortices related to various solutions of 
the equation above can in principle be only inferred by a direct study of equations of motion. 



It is possible to determine all the solutions of Eq. (3.13), for example using an explicit 
parameterization of the elements of SU(3). In what follows we will analyse three types of 
solutions that cannot be related by SU(2) color-flavor transformations: 1) "BDM" case, 2) 
"CP 1 " case, and 3) "Pure color" case: 

1) "BDM" case 

The first possibility is a closed loop generated by T 8 in SU (3) and the electromagnetic 
T EM alone: 

e il 8 (-fi)T 8 _ e -i8( v ) e -ia( v )T™ 

4 if = 2tt 

7 8 (2vr)/ v / 6 + a(27r)/3 = -2tt/3, #(2tt) = 2tt/3 . (3.14) 

The equation above determines the phases 7 s and a only up to a linear combination. This is 
a consequence of the fact that the two generators T 8 and T EM are proportional and indistin- 
guishable from each other on diagonal configurations. The configuration above is invariant 
under SU(2) color-flavor transformations. 

2) "CP 1 " case 

The second possibility is a closed loop generated by winding in the SU(3) group around 
the T 3 direction too, in addition to T 8 and T EM . 

e i(7 3 (v)T 3 + 7 8 ( ¥3 )T 8 ) _ e -ie( V ) e -ia&)T™ 

4 if = 2tt 

7 8 (2vr)/ v / 6 + a(27r)/3 = vr/3, 7 3 (2vr) = iv^vr, 9(2n) = tt/3. (3.15) 

The configuration above will not be preserved by color-flavor transformations, and we can 
generate a orientational moduli space using SU(2) color-flavor transformations. In fact, we 
have the most general configuration of this type in the form: 

7 3 (27r)T 3 ->• 7 6 (27r)T b |7 6 | = V2tt, 6=1,2,3, (3.16) 
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where — > denotes a replacement; the T a above are the generators of SU(3) that commute 
with T 8 and form an SU(2) subgroup. 
3) "Pure color" case 

In terms of the vector |7 a | introduced above, the two previous cases correspond to |7 a | = 
and |7 a | = V2tt. These are the only two cases for which a non-trivial electromagnetic phase 
is allowed. In all the other cases, we must have: 



e il a {ip)T a _ e -i9(<p). 

JJ. tp = 2tt 

#(2tt) = 2tt/3, a(2vr) = 0. (3.17) 

We called this case "pure color" because it doesn't involve electromagnetic transformations, 
and is thus equivalent to the well-known case of pure color vortices without the electromag- 
netic coupling. Notice that this case spans a whole CP 2 , while SU(2) color-flavor transfor- 
mations generates only CP 1 orbits (apart from the case where only 7 8 is non-zero, which is 
invariant). 

We stress again that the cases listed above are just a consequence of boundary conditions 
when we search for closed loops at spatial infinity. Moreover, all the cases are topologically 
equivalent. In the next subection, we will study how the cases above translate into stable 
vortex configurations. 



3.2.1 BDM vortex 



The "BDM case" corresponds to the vortex studied by Balachandran, Digal and Matsuura 
in Ref. (9j. Since we only need T 8 to generate the correct winding for this vortex, the idea 
is to restrict the action (3.2) to include only the gauge fields A 8 and A em , by keeping all 
the other gauge fields to zero. Formally the action reduces to that of a (7(1) x £7(1) gauge 
theory, which we can then express in terms of the massless and massive combinations (3.7): 

SU(3) C x E7(1) BM E7(l) 8 x U(1) EM ~ U(l) x U(1) M , (3.18) 

where the arrow means that we truncate the model to its Abelian sub algebra. The La- 
grangian reads 



C = Tr 



Tr 



1 : * /.•" pOy ... 

4 2 tJ 

__ F M F Mij + j^v^ty^ - A 2 ($ t $) 2 + m 2 &<& 
4 2 J 



A^Tr^]) 2 - 



3m 



4(3Ai + A 2 )' 
(3.19) 



with 



MrpM 



V i = d i -ig M A M T 



9m 



= ^e 2 + 3^ 2 /2, T M =T 



EM 



(3.20) 



which is applicable only for diagonal configurations for the order parameter $. In the action 
above, the massless field decouple completely, as expected, while the massive field Am couples 
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to $ in the standard way. The BDM vortex is constructed analogously to (2.14) with the 



only difference from the un-coupled case being the new coupling constant gu instead of g s : 



:,Kf 



AMrpM 



A r 



f(r) 








— ^-[l-h(r)]T M , 
9m r 2 



A° = 0. 



(3-21) 



As mentioned in Section 2, the tension of the vortex decreases monotonically with the gauge 
coupling. Since we have > g s , the BDM vortex has a smaller tension as compared to the 
corresponding vortex in the un-coupled case. 



3.2.2 CP 1 vortex 

The CP 1 case is considered for the first time in this paper. We have identified a vortex 
configuration which is a solution of the equations of motion that corresponds to these new 
boundary conditions. To see this, we notice that, similarly to the previous case, we can 
consistently restrict the action to include only the gauge fields corresponding to generators 
commuting with T 8 . This corresponds to formally reduce to the case 

SU(3) C x U(1) EM SU(2) C x Z7(l) 8 x U(1) EU ~ SU(2) C x U(l) x U(1) M , (3.22) 

similarly as what we have done for the BDM case. The truncated lagrangian now is the 
following: 



C = Tr 



Tr 



1 3 

4 2 ij 
3m 4 



+ 



4 J 



m 



Ai(Tr [$ t $]) 2 + 



4(3Ai + A 2 ) ' 



V i = d i -ig M A M T M -ig 8 A b I*, 6 = 1,2,3; [T 6 ,T 8 ]=0, (3.23) 

where the index b is relative to the SU(2)c factor. As before, the massless combination 
decouples completely, and we can set it to zero. The simplest vortex configuration of the 
type considered here has the following diagonal form: 

/ g x {r) 
*(r,<p) CP i = A CFL e^f(r) 
V g 2 (r) 

a m t m = _IJ_^ tl _ h {r)\T M , 
2 gu t 

A*T* = ±\ e -^[l-l{r)]T\ (3.24) 
y2g s r 
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As explained in Eq. (3.16), we can generate a full CP 1 of solutions by applying SU(2)c+f 
rotations to the configuration above. Once the ansatz (3.24) is inserted into the equations 
of motion (2.13), we get the equations written in the Appendix. We have numerically solved 
these equations, and determined the energy of the vortex configuration. The tension, as 
schematically shown in Fig. pi is higher for the CP 1 vortices than for the BDM vortex. 



This result can be intuitively understood if we recall the observation of the previous Section, 
where we noticed that the tension of a color vortex decreases monotonically with the gauge 
coupling. The CP 1 vortex is built with both g s and gu couplings, thus, the interactions 
depending on g s < gu contribute to increase the tension with respect to a vortex built 
exclusively from g^- 

Since the CP 1 vortex is a solution of the equations of motion, it is a critical configuration 
for the energy density functional. In Section [4], we will argue that it must be a local minimum, 
and thus a metastable configuration. 



3.2.3 Pure color vortex 

Let us now come to the third case, the "pure color" case. This situation is realized when 
the boundary conditions are the same as the case wihtout electromagnetic coupling. This 
means that at infinity, only the non-Abelian gauge fields gives non-trivial winding, and thus 
non-vanishing fluxes, while the electromagnetic gauge field is zero everywhere. However, this 
situation is generically not compatible with the coupled equations of motion. An easy way to 
see this is the fact that boundary conditions in the pure color case imply a non-zero flux for 
the massless field Since A® is massless and unbroken, there is no topology, compatible 
with the equations of motion, stabilizing and confining the flux. We thus expect that no 
solutions exist in general for the pure case. However, a special configuration that we call 
"pure color vortex" is an exception to this statement. This configuration is defined as the 
one corresponding to the following boundary conditions: 

e <7»(*r)T« = e ~i0(2n) . yj = Q (3 35) 

Notice that the conditions above are perfectly consistent. We can take any configuration of 
7 a phases and set 7 s to zero using a gauge (-flavor) transformation. This means that we can 
consistently set: 

A 8 i= A^ m = 0, Af = A° = 0. (3.26) 

In fact, A e ™ = means that for the pure color vortex we can consistently restrict the action 
by simply dropping all the terms involving Af m . As a consequence, the pure color vortex is 
exactly the same configuration as we get in the un-coupled case. Moreover, it satisfies the 
full equations of motion of the coupled case, because of the consistent restriction. Notice 
that the dependence of the electromagnetic gauge coupling disappears also completely from 
the restricted action. This means that the tension of the pure color vortex is also the same 
as the tension of the un-coupled vortices, involving only g s . As represented in Fig. [3j the 
tension of this vortex is larger than both the BDM and the CP 1 vortices. 

Because of this very fact we are lead to conclude that the pure vortex is in fact a stationary 
point for the energy functional, but it corresponds to a local maximum. Both the BDM and 
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BDM vortex ~ r 9 




Tcpi 
Tbdm 



pure color 
vortex 



Figure 3: The three different special vortex solution explained in the text. The BDM 
vortex sits in the bottom left corner of the toric diagram and has the lowest tension. 
The other two diagonal vortices are located at the remaining corners and are connected 
by co lor- flavor transformations. All the upper right edge of the diagram represents the 
whole family of vortex configurations related by this SU(2)c+f transformations, namely 
a CP 1 moduli space. This vortices have a slightly larger tension compared to the BDM 
configuration. The dashed line in the middle of the diagram corresponds to vortices not 
winding along T 8 direction in color group. These vortices have the largest tension and 
do not involve electromagnetic gauge field. We reported a qualitative behaviour of the 
potential between these solutions under the diagram. 
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the CP 1 vortices are on the other hand local minima, with the BDM vortex being the 
absolute minimum. The CP 1 vortices are however metastable, and if long-lived can play a 
crucial role to the physics of the CFL superconducting phase together with the BDM vortex. 
The whole situation is schematically summarized in Fig. [3} Numerical values of the tensions 
are compared as 

T pure - r BDM = 0M8MeV 2 , T pure - T CP i = 0.0066MeV 2 , (3.27) 

for the same "realistic" choice of parameters made in Fig. [2j where in addition we have chosen 
m G = 92MeV. 

We conclude this Section by recalling that the low energy physics of an uncoupled color 
vortex is described by a CP 2 non-linear sigma model, as shown in Ref. |15|. When we 
consider the effects of the electromagnetic coupling, the effective theory should become a 
U(l) gauged CP 2 model (X9j , where we have to include also the effective potential sketched 
in Fig. |3j 



3.3 Magnetic fluxes 

There are two main differences between color magnetic flux tubes with and without electro- 
magnetic coupling. The first one, as we have already examined, is the lifting of the moduli 
space CP 2 to leave the stable BDM vortex and the family of metastable degenerated CP 1 
vortices. The second is the fact that coupled vortices now carry a non-trivial electromagnetic 
flux. This is given by the fact that coupled vortices are made of the massive field Af' 1 , which 
is in turn a linear combination of color and electromagnetic gauge fields. 
The BDM vortex carries a quantized Af 1 flux: 

d>f DM = lA M -dl = — . (3.28) 
J 9m 

Because of the mixing in the ground state, this means the following non-quantized fluxes for 
the color and electromagnetic fields: 

r /2 1 2tt vm 5 2 2tt r9 2e 2 , 
<f> 8 — , / $ EM = A 2 = (3 29) 

BDM 1 o 1 i 19 ' BDM -i | ro ' — o 9 ' 

V 3 1 + o z g s 1 + o z e 3 gj 



The fluxes of the CP 1 vortices can be similarly determined: 

$ 3 Cpl+ = ^A 3 -dl = ^=^, $ 3 cpl = -^ CPi . (3.30) 

Moreover, the quantized circulations of the BDM and CP 1 vortices are the same as that of 
usual un-coupled vortices, and equal to c.b/3. 

Notice that the expressions determined above imply that a color-neutral bound state of 
vortices necessarily carries also no electromagnetic flux, and vice versa. In particular, in the 
un-coupled case we need at least a bound state of three vortices, (f, g, b) to obtain a color-less 
state, which is nothing but the U(1)b vortex. In the coupled case, this "minimal" color-less 
bound state is obtained with the combination (BDM, CP|, CPl). The bound state carries 
an integer circulation c# = 1. 
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r 9 b 

\r 4" "\r 

BDM CP+ CPl 

Figure 4: The three different types of vortices we consider. The BDM vortex is labelled as 
anti-red f vortex as explained in the text, and carries twice electromagnetic flux compared 
to the other solutions, the anti-blu b and anti-green g vortices. 



3.4 Quantum mechanical decay 

As we have seen, the CP 1 vortices are clasically metastable because of the potential barrier 
by the presence of the pure color vortex. However here we show that the CP 1 vortices 
are quantum mechanically unstable, and they can decay to the BDM vortex by quantum 
tunneling. We estimate the decay probability. 

As discussed in the previous subsections, we do not expect static solutions of the equations 
of motion apart from the BDM, the CP 1 and the pure color vortices. However, we can try 
to define an effective potential interpolating between the three types of solutions. At generic 
orientations, boundary conditions are the same as those for the un-coupled vortices. In 



fact, exactly the un-coupled vortex evaluated on the coupled action Eq. (3.2) gives the same 
tension as that of un-coupled vortices. This value of the tension is an upper bound for 
vortex configurations that have a fixed boundary condition corresponding to a generic point 
in CP 2 . The tension of the configuration that really minimizes the energy, for that fixed 
boundary conditions, defines an "effective" potential on the CP 2 moduli space induced by 
the electromagnetic interactions. Moreover, since the tension of vortices is mainly modified 
by the contribution of the mixed Af 1 , which couples with a larger gauge coupling, we expect 
the qualitative behavior of the effective potential to be of the form represented in Fig. [3] 

This qualitative picture is enough to address some important features of coupled color- 
magnetic flux tubes. The fact that the potential has more than one local minimum allows 
for the existence of kinks interpolating between the two vortices corresponding to the various 
(meta) stable configurations; these kinks are interpreted as confined monopoles. Moreover, 
the presence of kinks, and the higher tension of CP 1 vortices with respect to the BDM vortex, 
implies a decay rate of the former vortices into the latter through quantum tunneling. This 
tunneling proceed by enucleation of kink/ anti-kinks pairs along the vortex. These problems 



are very similar to those already studied for pure color vortices in Ref. 17 , where the 
potential was generated by quantum non-perturbative effects of the CP 2 non-linear sigma 
model. 
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We have a typical scale of the potential given by 



— Tpure Tbdm > (3.31) 

that has been evaluated numerically (see Eq. (3.27)). We must have semiclassical kinks with 



masses of order: 

M kink ~^/AV. (3.32) 

We are now ready to estimate the decay probability following Ref. [26] . The enucleation 
of a couple of kinks costs an energy of order M k i nk ~ V ' AV . Moreover, they are created at a 
distance L such that the energy cost for the pair production is balanced by the energy gain 
due to the presence of a intermediate BDM vortex with less energy: LAV ~ VAV. The 
decay probability rate per unit length is thus: 

P „ e -LVAV „ e -l ^ (3> 3 3) 

which is a quantity of order 1. Notice that this estimate is fully justified only in the case 
the size of the kink is negligible as compared to the string size, while in our case the two 
sizes are comparable. This estimate is very rough, and has to be corrected including possible 
Coulomb interactions between the kink/ anti-kink pair. However, it should correctly capture 
the order of magnitude of the decay probability. It then turns out that CP 2 vortices decay 
almost instantaneously in the CFL phase. This implies that the CFL phase with included 
electromagnetic interactions can support only a particular type of color flux, the one we 
called f which is carried by the BDM vortex. 



4 The comparison with the other potential terms 

In this section, we compare the potential generated semi-classically by the electromagnetic 
interactions with the other potentials, i.e., the quatum mechanically induced potential and 
the potential induced by the strange quark mass. 



4.1 Quantum mechanical potential 



Our task here is to compare the qualitative potential generated semi-classically by the electro- 
magnetic interactions to the quantum potential generated by CP 2 . The quantum potential 
has three vacua given by Ref. |27| 



E CP 2(k) ~ 3Ac p2 



1 + const 



2-Kk 



k = 0,1,2, 



(4.1; 



where 9 is a quantum generated 9 angle. The A scale of the sigma model is given by 



m 



exp 



3 A, 



(4.2) 



The constant c has been evaluated numerically in Refs. [17] and it is found to be of order 1. 
The vacua above are considered as coming from a quantum potential that oscillates between 
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the minima above and barriers of potential of height A<cp2. This implies the existence of 
kinks of mass: 

M kink ~ A cp2 . (4.3) 

The formula above has been evaluated directly in terms of the coefficients of the Landau- 
Ginzburg model. If we use the formulas derived at weak coupling, we also obtain K\j\\ ~ 
/x 2 /T c <C 1. The numerical value of Acp2, is thus exponentially small and negligible as 
compared to the semiclassical potential induced by electromagnetic interactions. 



4.2 Strange quark mass effects 

So far we have considered the CFL phase at very high densities, where quark masses can be 
neglected, including the mass of the strange quark. The effect of a non-zero mass for the 
strange quark corresponds to the introduction of the following effective potential term in the 



Landau-Ginzburg description Eq. (3.2) |28|: 

\'„, = r'IY <!> T ( jjl 2 + T 6 ) <!> 



m 2 \i 



The effects of this potential on the CP 2 moduli space of pure color vortices have been 
studied in Ref. |16|. It can be considered as generating the following effective potential 

V CP 2 = e J dx 2 Tr [$tT 3 $] = Dd^p _ 1^2). (4. 5 ) 
Where D is determined in terms of the vortex profile functions: 

/■oo 

D = vreA e / drr {f - g 2 ) , A e = A CFL (m 2 m 2 + 2e/3) , (4.6) 
Jo 

and |0i| 2 + |02| 2 + |03| 3 = 1 parameterize the CP 1 s which, in the toric representation of CP 2 
are all parallel to the long side of the triangle. In our notations, f = (1,0,0), g = (0, 1,0) 
and b = (0,0, 1). The vortex with less energy turns out to be then the one with |0 2 | — 1; 
or the g = (0, 1,0) vortex. When we turn on the electromagnetic interactions, this vortex 
corresponds to what we called the CP| vortex. We then see that the introduction of strange 
quark mass makes the BDM vortex unstable through decay to the CP| vortex. Our new 
solution is thus the one that can play the main role in the CFL phase at intermediate 
densities, when the mass of the strange quark cannot be neglected. Moreover, it removes the 
degeneration between the CP 1 vortices. We have compared the two effects with numerical 
simulations, and the results are: 

T = 5897MeK 2 , D = 4231Me^ 2 , AV = 0.0275MeV, (4.7) 

with the usual choice of parameters and in addition m s = 150MeV. 

We see that the potential induced by electromagnetic interactions is comparable with 
the effects of the strange quark mass only when the strange quark mass is very negligible as 
compared to the chemical potential (e m 2 ). These densities are not realistic in neutron 
stars core, for example. We then have found that the CPf, with its different flux with respect 
to the previously studied BDM vortex, is the the most relevant to the physics of neutron 
stars in their inner core. 
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5 Conclusions and Discussion 



The non-Abelian vortex in the CFL phase has the degeneracy of the Nambu-Goldstone modes 
CP 2 , associated with the color-flavor locked symmetry SU (3) c+f spontaneously broken down 
to a subgroup in the vicinity of the vortex, when the electromagnetic coupling is neglected. 
Once the electromagnetic coupling is taken into account, the flavor symmetry SU(3)f is 
explicitly broken. We have studied the electromagnetic coupling of the non-Abelian vortices 
and have found that the degeneracy of the CP 2 modes is removed with the effective potential 
induced. We have found that the stable BDM vortex solution and a family of metastable vor- 
tex solutions parametrized by CP 1 , both of which carry ordinary magnetic fluxes in addition 
to color magnetic fluxes. We have discussed quantum mechanical decay of the metastable 
CP 1 vortices by quantum tunneling. We also have compared the effective potential which 
we found with the other known potentials: the quantum mechanically induced potential and 
the potential induced by the strange quark mass. We have found that it is more dominant 
than the former but is negligible compared with the latter. Therefore all vortices including 
the BDM vortex decay into one of CP 1 vortices. 

The results we have found are important for the study of the physics of neutron stars. The 
very high rotation speed of these compact stars causes the formation of the non-Abelian semi- 
superfluid vortices we analyzed. These solitons arrange in a lattice and they are constrained 
to be together to cancel the total color flux. Each vortex carries also a non-zero magnetic 
flux due to electromagnetic coupling, but as summarized in Fig.|4]the total flux of this vortex 
bound state is zero. Then we can say that no magnetic field generated by vortex formation 
comes out from the CFL core. However, it is well known that neutron stars possess a huge 
magnetic field, which should penetrate in the core as well. Thus there should be some 
non-trivial interaction between the external magnetic field spreading into the core and the 
semi-superfluid vortices present there and it is thus an interesting problem to know how the 
system behaves in this kind of situation. 

Problems similar to those considered in this paper can be analyzed in supersymmetric 
models with non-Abelian gauge group 10 , IT) , where a global color-flavor locked symmetry 



is unbroken in the vacuum. Interesting physics can be found when gauging the whole or a 
part of the color-flavor group |29| 
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2 X Cj/jX^ 



$(r,y?) f = A CFL | g(r) | , AfT 8 = J -- [1 - h(r)\ T 8 , (A.3) 



A Appendix 

Phases and Fluxes 

A generic boundary condition for the order parameter, representing a closed loop in the 
global-gauge groups can be represented by the following phases 

($(oo, 2tt)> = e ld e^ aTa ($(oo, 0)) e iaT ™ , (A.l) 

which correspond to the following fluxes 

i 2tt g s r 2 ' i 2ie r 2 

$ a = 2^ $™ = 2 a! i (A. 2) 

Equations of motion for the un-gauged and BDM vortices 

With the usual ansatz 

e^/(r) 
9(r) 

g(r) J V 69s r 

we obtain the usual equations for the profile functions f(r), g(r) and h(r): 

f" + f i ~ ®££f ~ ff(f + 2g> - 3) - ff(f - gi) = 0; 
g" + 9 i- { -^g - r -fg(P + 2g 2 - 3) + fg(P - g 2 ) = 0; 

h" - £ - f (<? 2 (/i - 1) + / 2 (2/> + 1)) , (A.4) 

with 

2 o 2 a 2 T y 2 ^/i 2 2 4A2A 2 , 2 n /a r\ 

™ 5 = gs cfl 1 ) m * = m x = — ^ — ' m ^ = • ( A - 5 ) 

The equations for the BDM vortex are obtained by just substituting T 8 — > a/2/3T m and 
BDM vortex in terms of A? and 

The vortex with minimum energy, including the BDM vortex, is constructed in a background 
where the massless gauge field A® vanishes. We then have for the corresponding phases: 



-\j¥t a + fj 8 = 0; 7 8 /V6 + a/3 = 2tt/3; 

7 - y 3 i +5 2 a — 1+J2 , _ 3 g 2 . \ / 



Then we have 
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Equations of motion for the CP 1 vortices 

By substituting the following ansatz: 



CP± 



9i(r) 
A,,, i e^f(r) 





MrpM 



ATT 



1 1 eijX 3 
2g M r 2 



[1 - h(r)]T M , A\T 



i 3t-i3 



1 1 f--r 3 

13 \l-l{r)]T 



V29s r2 



into the equations of motion (|2.13|) we obtain: 

7 



r + f 

9'{ + 9 i 

g'l+ B f- il+h/ t* /2l)2 92 - fgi(P + 9l + 9l - 3) 



(1+V2+3/2Z) 2 
9r 2 



(fe-i) 2 „ _ !^ 
9r 2 yi 6 



h" - 



I" 



7/(/ 2 + 3? + 9l - 3) - ^/(/ 2 - gf/2 - ^ 2 /2) = 
Si(/ 2 + 9l + S 2 2 - 3) - f gi (g 2 - f 2 /2 - gl/2) = ( 

f 9 2{gl-f 2 /2-gl/2) 
(2gl(h - 1) + / 2 (1 + ty2 + 3/2/) + + V2 - 3/2/)) ; 
I - f (/ 2 (1 + ty2 + 3/2/) - g\{\ + h/2 - 3/2/)) . 
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